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Abstract 

A phenomenological model for the nucleon structure functions is presented. 
Visualising the nucleon as a cavity filled with parton gas in thermal equi- 
librium and parametrizing the effects due to the finiteness of the nucleon 
volume, we obtain a good fit to the data on the unpolarized nucleon structure 
functions. 
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Recent experiments have revealed some remarkable features of the nucleon structure 
functions Ff'". Data on deep inelastic scattering of muons off proton and deuteron targets 
|jlj] show that the quark sea in the nucleon is not flavor- symmetric, u(x) ^ d(x); the Gottfried 
sum |S g eJ (F 2 p - F 2 n ) (dx/x), at Q 2 = 4 GeV 2 , has the value 0.235 ± 0.026 compared 
to the usual quark model prediction of 1/3. This result has been confirmed by the observed 
asymmetry in Drell-Yan production of dileptons in pp and pn collisions ||. Most notably, 
the HERA electron-proton scattering data @] reveal a rapid rise of the proton structure 
function F%{x) as x decreases. Indeed over a wide range of small x, data from the various 
groups HH, for fixed Q 2 , are all well described by a single inverse power of x. Figure 1 is a 
log-log plot of the data on F%{x)/x (the combination that enters Sq) versus x. We see that, 
for fixed Q 2 , the data fall on straight lines defined by 

*IM = _£_ , (0.0004 < x < 0.2). (1) 

For instance, at Q 2 = 15 GeV 2 , the best-fit parameters are c = 0.229 ± 0.005 and m = 
1.22 ±0.01 §. 

Global fits to the nucleon structure data involve parametrizing the various parton densi- 
ties at some low Q 2 and evolving them to higher values of Q 2 relevant to observations. The 
fits so obtained |7j] have very high precision but contain several (typically ~15-20) arbitrary 
parameters and provide little physical insight into the structure of the nucleon. On the 
other hand, phenomenological models could give us some valuable clues into the physics of 
parton distributions in the nucleon. From this point of view the statistical models of the 
nucleon structure functions || have been quite interesting due to their intuitive appeal and 
simplicity. 

We present here a phenomenological model for the unpolarized nucleon structure func- 
tions by regarding the nucleonic contents as constituting a gas of noninteracting partons in 
thermal equilibrium. An attractive feature of this general framework is the natural expla- 
nation of the violation of the Gottfried sum rule: the excess of w-quarks over d-quarks in 
the proton implies unequal chemical potentials, hence unequal uu and dd seas, which leads 
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to So ^ 1/3. However, the ensuing structure function ^(x) vanishes like — > 0, in 
violent conflict with the data. To remedy this we invoke corrections arising from the finite- 
ness of the nucleon volume, by multiplying the parton density of states by a factor having 
inverse powers of the radial dimension, [1 + 0(1/ R 5 )]. We find that the small-x rise and 
other features of the nucleon structure functions are reproduced quite well. 

The Model 

We picture the nucleon (mass M) to consist of a gas of massless partons (quarks, anti- 
quarks and gluons) in thermal equilibrium at temperature T in a spherical volume V with 
radius R. We consider two frames, the proton rest frame and the infinite-momentum frame 
(IMF) moving with velocity — —1) along the common z axis. Our interest lies in the 
limit when the Lorentz factor j = (1 — v 2 )^ 1 ^ 2 — > oo. The invariant parton number density 
in phase space || is given by (quantities in the IMF are denoted by the index i) 



dn l dn g 



d 3 p % d 3 r l d 3 p d 3 r (2ir] 
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1 



exp[[3(E - a*)] ± 1 



= f(E), (2) 

where (3 = T _1 , g is the degeneracy (g = 16 for gluons and g = 6 for q or q of a given 
flavor), (E,p) is the parton four-momentum in the proton rest frame and f(E) is the usual 
distribution for noninteracting fermions or bosons. In terms of the Bjorken scaling variable 
x = pl/(Mvj), the phase space element can be expressed as 

d 3 p i dV= 27rpUpUMvydx) — 
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Ev 

2ti[Mxv 3 + -^-}dE Mdx d 3 r, 

sy2 



where for 7 — > 00 the expression in square brackets becomes Mx. For fixed x the parton 
energy E varies between the kinematic limits Mx/2 < E < M/2, where the lower limit is 
attained when p l T = 0. Consequently the parton number distribution dn % jdx in the IMF is 
simply proportional to an integral of the rest-frame distribution f(E): 



r M/2 

dn l /dx = 2WM 2 x dEf(E), (3) 

JxM/2 



where the factor V results from d 3 r integration. The structure function F 2 (x) is given by 



dril dn% 
dx dx 



F 2 (x) = xJ2< 
q 

The number distribution vanishes linearly as x — > (and also as x — > 1) and leads to 
the behavior of the structure function F 2 {x) ~ % 2 at small x, which disagrees with the 
observations noted in Eq. (1). 

In order to obtain the rise of F 2 (x) at small x, we shall modify the model to reflect effects 
arising from the finiteness of the nucleon volume V. Various studies of finite-size corrections 
(FSC) show that they are sensitive to the precise shape and size of the enclosure, the type of 
boundary conditions imposed on the wave function, and to the details such as whether the 



particles are strictly massless, whether chemical potentials are nonzero, etc. ||10|| . Moreover, 
these studies invariably involve some simplifying assumptions and thus their use is difficult 
to justify in the present context. 

In keeping with the phenomenological nature of the model, we have chosen to parametrize 
the correction due to the finiteness of the nucleon volume. This is implemented through the 
use of the dimensionless combination 1/(ER). We have chosen two alternative forms of 
parametrization, a form prompted by the empirical observation in Eq. (1): 

*' = 1 + (W' (4) 

and a general power series expansion in the variable 1/(ER): 

* 2 = 1 + is + (i^ + -- (5) 

where B, S(> 0), a, b, ... are arbitrary constants. We multiply the integrand in Eq. (3) by 
the function $ (= <£>! or $ 2 ) m order to incorporate the finite-volume effects in our model. 

The model described above is assumed to hold at a certain input momentum scale Q$, 
and if necessary can be evolved to higher Q 2 by means of the standard techniques in quan- 
tum chromodynamics (QCD). To complete the statement of the model, we demand the 
thermal parton distributions to obey the following three constraints at the input scale. The 
constraints on the net quark numbers in the proton are n u — n a = 2 and — = 1, i.e., 



VM 2 f 1 rM/2 f 6 6 1 

dxx dE — — -— $(ER) 



(2ir) 2 Jo JxM/2 [exp \J3(E — fi a )] + 1 exp [f3(E + p a )} + 1 

= n a — n & - (a — u,d) (6) 

Obviously, chemical potentials for heavy flavors are necessarily zero. As regards the third 

constraint, we assume that the longitudinal momentum fractions in the u, d flavors and the 

gluons add up to unity: 

VM 2 /-i , „ fM/2 t 6 6 

/ dx x I dE < 1 h 

(27r) 2 io Am/2 \exp[(3(E - p u )} + l exp [(3(E + /i u )] + 1 



exp [(3(E - fi d )] + 1 exp \J3(E + fi d )\ + 1 exp((3E) - 1 . 
The quark flavors s, c, ... which are not introduced in Eq. (7) show up at higher Q 2 as a 
result of QCD evolution. 

By interchanging the order of x and E integrations in Eqs. (6-7) and performing the 
^-integration analytically, we see that in order to keep the integrals finite, large powers of 
1/E are not allowed in the integrand. This requires that while using $x the exponent should 
be bounded, 5 < 3, and while using $ 2 only the first three terms can be present. Thus the 
model effectively has only two free parameters. 

To determine fj, u , fid and T, we solved the three coupled nonlinear equations (6-7) by 



the Davidenko-Broyden method Ull] . The resulting values of fj, u , fj, d and T are such that 
the left and right hand sides of these equations agree with each other to typically one part 
in 10 6 . The parton densities were evolved by means of the Gribov-Lipatov-Altarelli-Parisi 
equations Jjjf in leading order (LO), taking the input scale Ql = M 2 and Aq C d = 0.3 



GeV. Finally, the root-mean-square (rms) radius of the parton distribution was taken to be 



the same as the charge rms radius (p) of the proton; since p ~ 0.862 fm | 13fl , this yields 



R = .^573 p = 1.11 fm. 

Results and Discussion 

Since the two arbitrary constants B and S, or a and b in Eq. (4) or (5) are not known, we 
have determined them by fitting the deep inelastic scattering data on F%{x) at Q 2 = 15 GeV 2 
]. The results of our fit incorporating the finite-size corrections and QCD evolution are 



shown by the solid curve in Fig. 2. (Results presented here are based on $2; the alternative 
form $! gives an equally good fit.) Also shown for comparison in Fig. 2 are: (a) the (dot- 
dashed) curve labeled 'GAS' giving the prediction of the parton gas model which has no free 
parameters by virtue of the constraints, (b) the (dashed) curve labeled 'QCD' showing the 
effect of QCD evolution on the gas model, and (c) the (dotted) curve labeled 'FSC showing 
a fit to the data when only the finite-size corrections are introduced in the gas model. If $1 
is used in order to incorporate FSC, the fitted values of the two parameters are 

B = 0.269 and 5 = 2.14, 

and the corresponding temperature and chemical potentials are T = 63 MeV, \x u = 124 MeV 
and fid = 64 MeV. If, on the other hand, $ 2 is used, the fitted values of the two parameters 
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a = -1.88 and b = 2.24, 

and the corresponding temperature and chemical potentials are T = 72 MeV, \x u = 162 MeV 
and fi d = 81 MeV. 

To comment on the relative importance of the inputs, we focus on the curves in Fig. 2 at, 
say, x = 10~ 3 : a fit with FSC gives a very small value of F% ~ 0.02, reflecting the restrictive 
nature of the constraints. Leading-order QCD evolution does result in a value of F 2 which is 
significantly large but not large enough, F$ ~ 0.23. However, when the effects due to both 
FSC and QCD are included in the model, we obtain F$ ~ 1.1, which is consistent with the 
data. The presence of inverse powers of (ER) in $ is thus partially responsible for increase 
in F% at small x. 

As a test of the model, we show in Fig. 3 the prediction (solid curve) for the difference 
[F$ (x) — F^ipc)]. Also shown for comparison is the result (dashed curve) based on the 
parametrization of Gliick et al. M . The agreement with the NMC data is reasonable. 
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As for the Gottfried sum Sq, we have 



1 2 r 1 - 

3 6 JO 



1 2 VM 2 r 1 r M / 2 ( 6 6 

dxx dE\ -— - -— - \$(ER). 



3 3 (2tt) 2 Jo JxM/2 {exp\J3(E + n d )] + l exp \j3{E + fi u )\ + 1 

(8) 

The inequality Sg < | is thus a result of having in the proton, more valence u quarks than 
valence d quarks, {n u — riu) > (n^ — nj), implying that fi u > fid and hence the integral in Eq. 
(8) is positive. Our model predicts at Q 2 = 4 GeV 2 , the value So = 0.22 which is consistent 
with the experimental value Sg = 0.235 ± 0.026. 

The rapidity dependence of the W charge asymmetry in the reactions pp —> + ... is 



now known to a very high precision [flq] . It is a sensitive function of the quark flavor ratio 
d(x)/u(x) in the proton, in the range 0.007 < x < 0.24 at Q 2 = M^. The ratio u(x)/d(x) 
at < x >= 0.18 has been deduced to be about 0.51 by the NA51 collaboration ||. These 
and other predictions of the model, on the ratio (F^x) / F^x)), the quark and antiquark 
distributions q(x), q(x), q v (x) = q(x) — q(x) for various flavors, the gluon distribution g(x), 
the longitudinal momentum fraction carried by the charged partons, etc. will be given 
elsewhere |TB| . 

Now we briefly describe the salient features of some of the recent calculations of the 
nucleon structure functions, which use ideas from statistical mechanics. Mac and Ugaz [8a] 
calculated first-order QCD corrections to the statistical distributions and obtained a crude 
but reasonable agreement with F%(x) data for x ^ 0.2. The momentum constraint was not 
imposed and the fitted value of the proton radius (R) was 2.6 fm. Cleymans et al. [8b] 
used the framework of the finite temperature quantum field theory. They considered 0{a s ) 
corrections to the statistical distributions and obtained a good fit to the F^{x) data for 
x > 0.25. They also calculated the ratio <J L /ax in this region; it was a factor of 6 above 
the experimental value. Bourrely et al. [8c] considered polarized as well as unpolarized 
structure functions and presented a statistical parametrization (with eight parameters) of 
parton distributions in the IMF. Their framework allowed chemical potential for quarks as 



well as for gluons. The number constraints were not satisfied very accurately. QCD effects 
were not considered. xq(x) vanished as x — > and so it was not possible to reproduce the 
fast increase of the antiquark distributions for a; < 0.1. Bourrely and Soffer's [8d] approach 
was similar to that in [8c]. By incorporating QCD evolution of parton distributions and 
allowing the antiquark chemical potential to depend on x, they were able to reproduce the 
HERA data on Ff . 

In conclusion, it is noteworthy that the application of ideas of statistical mechanics to 
the point constituents of the nucleon can provide a simple description of all the observed 
features of the (unpolarized) nucleon structure functions down to the lowest x values so far 
explored. The model presented here has two free parameters which arise from our treatment 
of the finite-size corrections. 
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FIGURES 

FIG. 1. Log-log plot of the proton structure function data. Experimental data are from Refs. 
[4,5]; the error bars show statistical and systematic errors combined in quadrature. The straight 
lines are our fits described in Eq. (1), and are labeled by Q 2 = 15, 35, and 120 GeV 2 . Numbers 
have been scaled by the factors shown in parentheses for convenience in plotting. 

FIG. 2. Proton structure function F%(x) at Q 2 = 15 GeV 2 . Data points are as in Fig. 1. 
Solid curve is our best fit to the data. Also shown for comparison are: the (dot-dashed) curve 
labeled 'GAS' giving the gas model prediction, the (dashed) curve labeled 'QCD' showing the 
QCD-evolved gas model, and the (dotted) curve labeled 'FSC which is a fit to the data when 
finite-size corrections are included in the gas model (without QCD). 

FIG. 3. Difference (F^ — Fg) versus x, at Q 2 = 4 GeV 2 . Experimental data are from Ref. [1]; 
errors are statistical only. Solid curve is the prediction of our model. Dashed curve is based on the 
parametrization of Gliick et al. [7]. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9506367v2 
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This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9506367v2 



Figure 2 
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This figure "figl-3.png" is available in "png" format from: 
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Figure 3 




